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Abstract The present paper reconsiders the Newtonian limit of models of 
modified gravity including higher order terms in the scalar curvature in the 
gravitational action. This was studied using the Palatini variational principle 
in [Meng X. and Wang P.: Gen. Rel. Grav. 36, 1947 (2004)] and [Domi'nguez 
A. E. and Barraco D. E.: Phys. Rev. D 70, 043505 (2004)] with contradict- 
ing results. Here a different approach is used, and problems in the previous 
attempts are pointed out. It is shown that models with negative powers of 
the scalar curvature, like the ones used to explain the present accelerated 
expansion, as well as their generalization which include positive powers, can 
give the correct Newtonian limit, as long as the coefficients of these pow- 
ers are reasonably small. Some consequences of the performed analysis seem 
to raise doubts for the way the Newtonian limit was derived in the purely 
metric approach of fourth order gravity [Dick R.: Gen. Rel. Grav. 36, 217 
(2004)]. Finally, we comment on a recent paper [Olmo G. J.: Phys. Rev. D 
72, 083505 (2005)] in which the problem of the Newtonian limit of both the 
purely metric and the Palatini formalism is discussed, using the equivalent 
Brans-Dicke theory, and with which our results partly disagree. 

Keywords Newtonian limit • f{R) gravity • Palatini variation • metric 
variation 



1 Introduction 

The cosmological puzzle of the late time accelerated expansion of the universe 
has motivated a lot of authors to consider models of modified 
gravity for its theoretical explanation. One type of models is the so-called 
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fourth order gravity, where the gravitational action involves higher order 
terms in the scalar curvature, R (see for example 013 )• This idea is not new; 
it has been studied in the past for explaining early time inflation The 
main difference between the two approaches is that in order to generate early 
time inflation one needs positive powers of the scalar curvature in the action, 
whereas to explain the present accelerated expansion, negative powers of R 
should be present. 

In 10 VoUick considered a further modification, by using the action of [21 
together with the Palatini variational principle. According to this principle, 
the metric and the affine connections are considered as geometrically inde- 
pendent quantities. The action has to be varied with respect to both of them 
to give the field equations. This approach gives identical results to those ob- 
tained with the metric approach, if the standard Einstein-Hilbert action is 
used, the variation with respect to the metric giving the Einstein equations 
and the variation with respect to the connections giving the expressions for 
the Christoffel symbols However, differences arise between the results 
obtained with the two approaches when a more general action is used. 

The Palatini approach can be viewed as a further generalization of fourth 
order gravity. However, standard general relativity does not need the assump- 
tion that the connections are the Christoffel symbols of the metric; this can 
simply be derived from the variation of the action. In that sense the Palatini 
approach seems more appealing since it requires one less assumption than 
the metric approach. Another important aspect is that a gravitational action 
with higher order terms in R, treated in this formalism, gives a set of sec- 
ond order pde's plus an equation involving the connection, which is trivial 
to solve and interpret using a conformal transformation .10, . Remarkably, in 
vacuum, the theory reduces to standard general relativity plus a cosmologi- 
cal constant. Fourth order gravity in the metric formalism, however, leads to 
fourth order differential equations (hence the name), which are very difficult 
to treat. Thus, also from the practical point of view, the Palatini approach 
seems preferable. It has also been shown that stringent constraints can be put 
on models of cosmological interest in the metric formalism, like the one dis- 
cussed in |2j , via confrontation with the solar system experiments |12| . Even 
though this problem can be overcome by more sophisticated models this 
requires significant fine tuning of the various parameters. Additionally, Dol- 
gov and Kawasaki have shown in J13' that a violent instability appears in 
models like j2j in a weak gravitational regime within matter, when treated 
within the framework of the metric formalism, but the Palatini version of the 
same models is free of such instabilities (see |14|'). 

Theories like those presented in [Jj or ^01 seem to explain the present 
accelerated expansion |7lll4| but there are obscure points about their behavior 
in a weak gravity regime. Having the correct behavior in the Newtonian 
limit is, of course, fundamental for any proposed theory. In |15| constraints 
were obtained for the models in the metric formalism, in order for them to 
reproduce gravity as we know it from experiments in weak field. The picture 
in the Palatini formalism is less clear. In JS] Meng and Wang claim that 
all models with inverse powers of the scalar curvature in the action give a 
correct Newtonian limit. On the other hand, in ^3 it is claimed that this 
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is not true and that there are constraints on the form of the Lagrangian. 
During the preparation of this paper the post-Newtonian expansion of the 
discussed theories was studied in |18II19| both for the metric and the Palatini 
formalism, by using the equivalent Brans-Dicke theory. The results relevant 
to the Newtonian limit seem to be in agreement with those of \15\ and |17j . 
The line of reasoning in is basically no different from that used in JH| ; 
at least for what concerns this study. Thus, from now on, any reference to 
[TH] will actually include also |19| . 

In the present paper, we reconsider the weak field regime in the Palatini 
formalism and show that any reasonable model gives the correct Newtonian 
limit, following a different approach from both ^B], ^7] and jJHj. This is 
done in section |21 Section |3| is dedicated to detailing the points at which 
our analysis disagrees with that in JS| ■ The discussion in section El raises 
doubts about whether the studies of the Newtonian limit within the metric 
formalism performed by Dick ^S] and Olmo are viable. We comment 
on this in section^ Section |S1 contains conclusions. 



2 Palatini formalism 

We will consider theories of gravity coming from an action of the form 

^^hj '^^^^af{R) + Sm, (1) 

where f{R) is a polynomial including both positive and negative powers of 
R, K = SttGn, and Sm represents the matter action. The energy momentum 
tensor is given by 

In the Palatini formalism, when we vary with respect to the metric we get 
f'{R)R^, - = (3) 

and contracting, 

fiR)R-2fiR) = KT, (4) 

where the prime denotes differentiation with respect to R and T is the trace 
of the stress energy tensor. 

In and |17j the authors expand around de Sitter in order to derive 
the Newtonian limit. We write 

R = Ra + Ru (5) 

where Rq is the scalar curvature of the background de Sitter spacetime and 
Ri is the correction to Rq, including all possible terms, with Ri/Ro being 
considered as a small quantity. We will need to calculate /(i?o + Ri) and 
/'(i?o + The usual approached is to Taylor expand around R = Rq and 
keep only the leading order terms in i?i but we will show that this cannot 
be done in the present context because Ri/Rq is not small. 



4 



Thomas P. Sotiriou 



Take as an example the CDTT model j^, studied by Vollick in the 
Palatini formalism. Then 

f{R) = R-^, (6) 
and a ~ 10^^''(eV)^ ^ 10^'^'^m^^. Expanding we get 

f{R) = /(i?o) + f'{Ro)Ri + \f"{Ro)Rl + ■■■ (7) 
and using Q we get 

/(i?) = /(i?o) + (l + ^) i?i - \%^Rl + . . . (8) 

where now i?o = o- It is easy to see then that the second term on the right 
hand side of the above equation is of the order of Ri , whereas the third term 
is of the order of R\/a. Therefore, in order to truncate the third term one 
needs Ri ^ Rf/a or 

a>i?i. (9) 

Note that this is not any exceptional constraint. Rq = a and so this is the 
usual condition for being able to truncate non linear terms in a Taylor ex- 
pansion. 

Let us now return to eq. For the CDTT model (eq. ®) this gives 

i? = i (-kT ± v/k2T2 + 12a2) . (10) 

When discussing whether a theory has a good Newtonian limit, we are in 
practice checking whether the field equations reduce to a high precision to 
the Poisson equation, under certain assumptions: energy densities should be 
small enough so that there are no strong gravity effects and velocities related 
to the motion of the matter should be negligible compared to the velocity 
of light. At the same time energy densities should be high enough so that 
the system under investigation can be considered gravitationally bound ^. 
As a typical example of a density satisfying the above criteria we can take 
the mean density of the solar system, p ~ 10~"'^^gr/cm^ ^. This estimate now 
implies that \a/KT\ ~ lO"^!^ where T ^ ~p (c is taken to be equal to 1). 
The "physical" branch of the solution given in eq. seems to be the one 
with the plus sign in front of the square root. In fact, given that T < 0, this 
branch ensures that the matter leads to a standard positive curvature in a 
strong gravity regime. Then 

R kT-So^/kT (11) 

^ For example, in the cosmological constant model one could consider, even on 
non-cosmological scales, densities low enough so that the correction coming from 
the cosmological constant dominates with respect to the matter density in the 
Poisson equation. This of course would not imply that this model does not have a 
correct Newtonian limit 

^ One could imagine that this would be the density of a dust cloud which, due 
to some instabilities, collapsed to form the sub-structures of the solar system (sun, 
planets, asteroids, etc.) in a fairly Newtonian way. 
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and Ri ~ —kT ~ np. Thus a/Ri ~ 10^^^ and it is now evident that condition 
(9) does not hold for the typical densities related to the Newtonian limit. 

Note that the situation does not improve even if we choose the "unphysi- 
cal" branch of eq. H1U|) which has a minus sign in front of the square root. In 
fact in this case i?i ~ a{'ia/ kT — 1/3) so the correction to the background 
curvature is of the order a and not much smaller than that as required in 
order to truncate the higher order terms in the expansion eq. 

In |17j , this fact was overlooked and only linear terms in i?i were kept in 
the expansion of f{R) and f'{R) around Rq. In even though they notice 
it in the final stages of their analysis and they use it to truncate some terms, 
the authors do not take it into account properly from the beginning, keeping 
again only first order terms (eq. (11) of |16j for example). 

An alternative way to attack the problem of the Newtonian limit is the 
following. We already know from relevant literature ^01 that the connections 
are the Christoffel symbols of the metric 

h^.u = /'(i?)5M- (12) 

For the CDTT model then, and if we define e ~ a? /R'^, eq. (jSJ takes the 
form 

(1 + e)R^, - i(l - e)i?.g^, - kT^,, (13) 

and 

V = (l + e)5M-- (14) 
Due to eq. 10} e depends only on T. Combining eqs. and ((T^ we get 

Rfif ~ ^Rhf_,„ + e (j^^hf_,„ + Ri_,J^ = kT^„. (15) 

Note that up to this point no approximation or truncation was used. We have 
merely expressed the left hand side of the structural equation of spacetime 
with respect to quantities depending only on the /i^j^ metric, which is con- 
formal to gfj,i,. However, using eq. l(Tn|l and (fTT|l we see that e ^ 10^^^ if we 
use consider the mean density of the solar system as before and even less for 
higher densities. So the two metrics are practically indistinguishable in such 
cases, due to eq. H14|l . Thus we can use the h metric to derive the Newtonian 
limit. If we assume that 

V='7M^ + 'iM^' (16) 

then the first two terms of eq. H15|l will give the standard Newtonian limit 
and the last two terms will give a negligible contribution, since they are 
suppressed by the e coefficient. A deviation of the order of 10^"*^ is far below 
the accuracy of any known experiment. In fact, one can consider densities 
several orders of magnitude lower and still get corrections which will be much 
below experimental accuracies. 

A critical point is that here we assumed that the metric is flat plus a small 
correction instead of de Sitter plus a small correction. Note, however, that 
we are not claiming that we are expanding around the background or any 
corresponding maximally symmetric spacetime. We are merely asking for the 
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matter to account for the deviation from flatness which is the basic concept 
related to the Newtonian hmit. In any case, de Sitter is essentially identical to 
Minkowski for the densities discussed, and the important corrections to the 
metric come from the local matter, not from considerations of the universe 
as a whole. 

According to the above the CDTT Lagrangian in the Palatini formalism 
gives a perfectly good Newtonian limit. The approach can be extended to 
more general Lagrangians and it is safe to assume that all f{R) theories in 
this formalism will give a correct Newtonian limit as long as the extra terms 
are suppressed by small enough coefficients. This happens because the form 
of eq. Q is such that in all such cases, R ^ np for densities relevant to the 
Newtonian limit. 



3 Recovering the Newtonian Limit through the Equivalent 
Brans Dicke theory 

The problem considered here has also recently been studied in I18J using the 
equivalence between f{R) theories of gravity and Brans-Dicke theories. It 
is shown there that the Palatini formalism of any f{R) theory is equivalent 
to a Brans-Dicke theory of gravity with Brans-Dicke parameter uj — —3/2 
(for the metric formalism see section 0)1. Keeping that in mind, one can 
generalize known results to obtain the post-Newtonian metric. Olmo used the 
fact that, for approximately static solutions, one can drop terms involving 
time derivatives to the lowest order, and the metric can be expanded about 
its Minkowskian value. The final relations giving the post-Newtonian limit 
of the field equations where 



where V is the potential of the scalar field </) and i^{T) = \og[(f>/(f>o]. The 
subscript in (po, and in any other quantity henceforth, denotes that it is 
evaluated at T = 0. Note at this point that the normalization by 4>q in 
this definition is not required. Olmo probably just added (subtracted) the 
constant log(0o) inside the brackets on the left hand side of eq. (|17|l (eq. 
(|18|l l using the fact that it remains unchanged. Thus we are not going to use 
it here and we are going to refer to f^{T) just as i^{T) = log[(j)]. The solutions 
of eqs. (fT7|) and lfTS|l are 



-v2 [hlo - niT)] 



2<j) 

' Kp+V{(j)) 

_ 20 



(17) 



(18) 




hlAt,x) 




(19) 



(20) 
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where Mq = <j)o J (Px'p{t,x')/(l). The effective Newton's constant G and the 
post-Newtonian parameter 7 are defined as 

where My = k~^0o / d'^x' [Vb/^o — Note that the appearing in 

these relations in J2| is what we denote here as just k. Even though we agree 
with the approach followed to derive eqs. (|19|) and H2U|I and on their validity, 
we disagree with the line of reasoning used by Olmo to argue that models 
with inverse powers of the scalar curvature do not have a good Newtonian 
limit. We will demonstrate this using again the CDTT model (eq. ©). 

As stated in different words in ^Hjj if we define the Newtonian mass as 
Mn = J (Px'p{t,x'), the requirement for a theory to have a good Newtonian 
limit is that GMq is equal to GnMj^, where N denotes Newtonian and 7 ~ 1 
up to very high precision. Additionally, the second term on the right hand 
side of both eq. (|19|l and eq. H2()|l should be negligible, since it acts as a 
term coming from a cosmological constant. fi{T) should also be small and 
have a negligible dependence on T. The above have to be true for the range 
of densities relevant to the Newtonian limit, as discussed before. Using the 
equation that define V and (j) found in J5| (see also [7|) one can easily show 
that 

'^-1 + ^. (23) 
y((/)) = 2aKy^(l>-l. (24) 

Additionally, for T = 0, i? = Via so, 

00 = 4/3, (25) 
Vo = 2aK/V3. (26) 

For the densities we consider we can use the parameter e defined above. Then 

2 

V{(f>) = 2k^ = 2KaVi, (27) 

and Mv '-^ a. It is very easy to see using eq. ((23J, (QUli and l(?7|l that 

G«-^, (28) 
87r0o 

7 « 1, (29) 

and (/) « 1 plus corrections of order a or smaller, which is well above the limit 
of any experiment. 
Additionally 

n{T) = log[0] = log + log [1 -I- aV^^T^] . (30) 
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Vq is of the order of a which is a perfectly acceptable value and f2{T) is 
negligible at the densities being considered and decreases even more when 
the density increases. Therefore, our previous result are valid and theories 
including inverse powers of the scalar curvature have a correct Newtonian 
limit in the Palatini formalism. 

Of course this result contradicts those reported in ^HI , even though the 
approach followed there seems to be satisfactory. The main reason for this 
problem seems to be the following. In TS" the fact that f2{T) should have a 
mild dependence on T is used to obtain a constraint for the dependence of (j) 
on T (eq. (26) of ^H])- Taking a number of steps this constraint is turned into 
a constraint for the functional from of f{R) (eq. (37) of 18 ) and from that a 
conclusion is derived about the its possible nonlinearity. We disagree with this 
line of thought. Such inequalities constrain merely the value of the relevant 
quantity at the point where it is evaluated and not its true functional form. 
One could probably use them to make some assumptions about the leading 
order term but not to exclude any terms of a different form, as long as they 
are negligible with respect to the leading order for the relevant values of R. 
This, for example, is the case for the CDTT model discussed above. Any 
constraint placed by the Newtonian limit has to hold over a certain range of 
relevant densities (and consequently curvatures), and at not for all densities 
as implied in |18j . 

It should be clarified, of course, that the approach followed in I18II10| is 
broad enough to investigate, apart from the Newtonian limit of the mod- 
els discussed, also their post-Newtonian behavior and the confrontation with 
the solar system experiments. Solar system tests (light deflection, Shapiro 
time delay, etc.) are much more sophisticatedthan the Newtonian limit of a 
theory and at the same time very different. They do not examine gravitation- 
ally bound systems, but are essentially vacuum tests, in which the presence of 
matter {e.g. solar winds) has to be taken into account as a correction. There- 
fore the relevant densities can be many orders of magnitude smaller than 
those associated with the Newtonian limit. The approach which we followed 
here only considered densities relevant for the latter and our objections to 
|18II19 | are confined to the discussion about the Newtonian limit. For discus- 
sions of constraints derived using solar system experiments see, apart from 

[THmni, also pniiTni^ . 



4 Metric formalism 

As already mentioned in the introduction, the Palatini and the metric for- 
malism give rise to different theories. Starting with the same action, we 
get the following field equations within the metric approach: 

f'{R)R^, - i/(i?)5p. - V^V,/'(i?) + g^,V^f'{R) = nT,,,, (31) 
where is the covariant derivative, and = V^V^. Contracting we get 



f'{R)R - 2f{R) + 3\/^f{R) = kT^, 



(32) 
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For the CDTT Lagrangian, we get from 




18a2 



i?4 



(33) 



The above equation relates the scalar curvature with the stress energy tensor 
in a dynamical way. This makes any analysis a lot more complicated than in 
the Palatini approach, where eq. Q was just an algebraic equation. 

In Dick expanded around the de Sitter solution to derive the New- 
tonian limit in fourth-order gravity. In order to do so, he used a Taylor ex- 
pansion for f{R) and f'{R), keeping only first order terms in the correction 
to the scalar curvature 6R. As we showed in the previous section, truncating 
higher order terms is not something trivial when models with terms inversely 
proportional to R are considered. The condition to be fulfilled is SR ^ a for 
the CDTT model, which is the simplest and most typical example. Note that 
this condition has to hold both in the metric and the Palatini formalism. We 
have already shown that this condition does not hold in the Palatini version 
of the theory by using eq. UlUfl . The metric formalism is, of course, totally 
different, and one would have to use eq. to make predictions on the 
value of SR. The dynamical nature of this equation requires a different treat- 
ment from the one used in the previous section. What is clear, however, even 
without going any further in the analysis is that it not at all safe to assume 
that in all cases SR ^ a. No proof that this condition holds is given in 
and therefore, one has to assume that higher order terms in the expansion of 
f{R) and f'{R) were actually truncated ad hoc, which leaves doubts about 
whether the analysis presented there really applies in such models. This indi- 
cates that, the question of whether fourth-order gravity models with inverse 
powers of the scalar curvature give good Newtonian limits, is still not an- 
swered, at least not rigorously. However, since, as already mentioned, the 
analysis of this problem seems to be completely different from the one we 
used for the Palatini approach, we intend to address it separately in future 
work. 

As already mentioned, in |18| the Newtonian limit of fourth order grav- 
ity was reconsidered using the equivalent Brans-Dicke theory. In the metric 
formalism, an /(i?) theory corresponds to an w = Brans-Dicke theory. In 
this approach no expansion of f{R) is formally needed. However, one has to 
expand the potential of the scalar field V{(j)) around (po to derive eq. (11) 
and (12) of ^Sl- is given in terms of f{R) by 



(see eq. (9) of JH])- Therefore, it is safe to assume that any problematic be- 
havior in the expansion of f{R) and f'{R) will be inherited by the expansion 
of V{4>). This indicates that our concerns, stated in the previous paragraph 
and involving the approach presented in |15j . remain relevant also for the 
approach followed in (l5) . 



v{4>) = i?/'(i?) - /(i?), 



(34) 
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5 Conclusions 

The Palatini version of modified models of gravity including higher order 
terms in the scalar curvature has been studied in a nearly Newtonian regime. 
This has been done mainly, using as an example the characteristic model of 
|1(J| . but the results obtained were shown to be general. It has been shown 
that such models can give a good Newtonian limit, as long as the coefficients 
of the higher order terms are reasonably small. Additionally, some of the 
present results seem to indicate that there is doubt on whether the approach 
used in |15| and involved first order expansions with respect to the scalar 
curvature, or equivalcntly the approach presented in 1181, are adequate to 
study the Newtonian limit of models with terms inversely proportional to the 
scalar curvature in the gravitational action, in the purely metric approach. 
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